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ABSTRACT 

We present preliminary results of a new quasar survey being undertaken with multi-colour optical imaging 
from the Canada-France-Hawaii Telescope. The current data consists of 3.83 sq. deg. of imaging in the i' 
and z 1 filters to a 10cr limit of z' < 23.35. Near-infrared photometry of 24 candidate 5.7 < z < 6.4 quasars 
confirms them all to be low mass stars including two T dwarfs and four or five L dwarfs. Photometric estimates 
of the spectral type of the two T dwarfs are T3 and T6. We use the lack of high-redshift quasars in this survey 
volume to constrain the z = 6 quasar luminosity function. For reasonable values of the break absolute magnitude 
M* 450 and faint-end slope a, we determine that the bright-end slope (3 > -3.2 at 95% confidence. We find that 
the comoving space-density of quasars brighter than M1450 = -23.5 declines by a factor > 25 from z = 2 to 
z = 6, mirroring the decline observed for high-luminosity quasars. We consider the contribution of the quasar 
population to the ionizing photon density at z = 6 and the implications for reionization. We show that the 
current constraints on the quasar population give an ionizing photon density <C 30% that of the star-forming 
galaxy population. We conclude that active galactic nuclei make a negligible contribution to the reionization 
of hydrogen at z ~ 6. 

Subject headings: cosmology: observation - galaxies: active - quasars: general 



1. INTRODUCTION 

After the combination of protons and electrons at redshift 
z ~ 1000 the hydrogen content of the universe remained 
largely neutral for some time. Observations of high-redshift 
quasars show that the neutral hydrogen in the diffuse inter- 
galactic medium (IGM) became fully reionized by z ~ 6 (Fan 
et al. 2004). In the standard picture (e.g. Barkana & Loeb 
2001) this reionization was the result of a rapid rise in the 
production of hydrogen-ionizing photons due to the birth of 
the first sizable population of stars and accreting black holes. 

A census of the star-forming galaxy population at z ~ 6 is 
now possible due to deep optical surveys carried out by the 
Hubble Space Telescope (Dickinson et al 2004; Bunker et al. 
2004; Bouwens et al. 2004). Whilst the luminosity function 
of these galaxies is now fairly well constrained (Bunker et 
al. 2004), translating this to the UV ionizing photon produc- 
tion rate relies upon a couple of uncertain factors: the spectral 
shape shortward of Ly a (which depends upon the stellar ini- 
tial mass function and metallicity) and the fraction of ionizing 
photons which escape the host galaxy. Therefore the ionizing 
photon production rate of the star-forming galaxy population 
remains quite uncertain and it is not clear whether this rate is 
high enough to reionize the universe at z = 6 (Stiavelli, Fall & 
Panagia 2004). 

An alternative source of ionizing photons is the quasar pop- 
ulation. Quasars have very hard spectra which makes them 
extremely efficient ionizing sources. It is possible that the 
earliest phases of reionization at z > 10 were carried out by 
low-mass accreting black holes also known as 'mini-quasars' 
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(Madau et al. 2004). The z = 6 quasar luminosity function 
is much more poorly constrained than the z = 6 galaxy lumi- 
nosity function and consequently the quasar ionizing photon 
output is very poorly known. The only quasars known at this 
redshift are from the bright Sloan Digital Sky Survey (Fan et 
al. 2004 and refs therein) and these luminous quasars make 
up only a tiny fraction of the total quasar luminosity density. 
The bulk of the luminosity density is emitted by quasars close 
to the 'break' in the luminosity function. 

Attempts to constrain the high-redshift quasar ionizing pho- 
ton output have therefore largely depended upon extrapolation 
of the SDSS quasar space density to lower luminosities us- 
ing parameters derived for the luminosity function at lower 
redshifts (Fan et al. 2001; Yan & Windhorst 2004; Meiksin 
2005). This necessarily involves large and uncertain extrap- 
olations. These studies show that, for the expected values of 
the luminosity function parameters, quasars fall short of pro- 
ducing enough ionizing photons for reionization. However, 
there is considerable debate as to how many photons are nec- 
essary for reionization since there is a large uncertainty in the 
dumpiness of the IGM which controls the importance of re- 
combinations (Madau, Haardt & Rees 1999; Haiman, Abel 
& Madau 2001; Oh & Haiman 2003; Iliev, Shapiro & Raga 
2005; Meiksin 2005). 

Dijkstra, Haiman & Loeb (2004) presented a novel method 
to constrain the quasar contribution to reionization by con- 
sidering the unresolved component of the soft X-ray back- 
ground. They found that, for their adopted quasar spectral 
shape between the UV and hard X-rays and dumpiness of the 
IGM, z = 6 quasars cannot provide enough photons to reion- 
ize the universe at 2a significance. Meiksin (2005) discussed 
this soft X-ray background constraint and showed that adjust- 
ing either the typical quasar spectral shape or clumping fac- 
tor within reasonable ranges produced a result consistent with 
reionization at z = 6 by quasars. 

Current constraints on the z = 6 quasar luminosity function 
also depend heavily on the properties of SDSS quasars. The 
luminosity distribution has been used to place a 2a range on 
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the bright-end slope of -4.2 < /3 < -2.2 (Fan et al. 2004). The 
lack of strong gravitational lensing amongst the SDSS sample 
can also place a constraint on (3 (e.g. Wyithe & Loeb 2002; 
Comerford, Haiman & Schaye 2002). Fan et al. (2003) and 
Richards et al. (2004) showed that this leads to a constraint in 
the range (3 > -3.5 to (3 > -4.5 at 2a significance depending 
upon the assumed values for the break magnitude and faint- 
end slope. Wyithe (2004) combined SDSS lensing constraints 
and the luminosity distribution of Fan et al. (2003) to find 
that (3 > -3.0 at 90% confidence. However, using the steeper 
slope preferred by the luminosity distribution of the expanded 
sample of Fan et al. (2004) would have the effect of weak- 
ening the conclusion of Wyithe (2004). A search for lower 
luminosity quasars at z > 5 in the Chandra Deep Field-North 
found one quasar at z = 5 .2 but none at higher redshift (Barger 
et al. 2003). 

The Canada-France High-z Quasar Survey (CFHQS) 
utilises the wide-area deep surveys now being undertaken at 
the 3.6m Canada-France-Hawaii Telescope with the 1 degree 
field-of-view camera MegaCam. In this paper we present the 
method for selecting z > 5.7 quasars and our initial results. 
The data used in this paper represent a very small fraction of 
the full area/depth that is planned. 

The paper is organised as follows. In Sec. 2 we discuss 
the expected optical/near-IR colours of point sources (quasars 
and stars) as observed with MegaCam and the near-IR cam- 
era CFHT-IR. In Sec. 3 we describe the optical and near-IR 
datasets and the search for high-redshift quasars. In Sec. 4 we 
discuss the limits on the high-redshift quasar luminosity func- 
tion which result from our work. Sec. 5 addresses whether 
z~b quasars can produce enough UV photons to reionize the 
universe. Cosmological parameters of Hq = 70kms _1 Mpc -1 , 
l^M = 0.3 and Ha = 0.7 are assumed throughout. 

2. SIMULATING THE COLOURS OF QUASARS AND STARS 

2.1. High-redshift quasars 

The sharp drop in flux across the Lya emission line in 
high redshift quasars makes colour selection the most efficient 
method for selecting quasars over large areas (e.g. Richards et 
al. 2002). This technique has been successfully utilised in the 
SDSS to discover all of the highest redshift quasars currently 
known (Fan et al. 2004). Spatially unresolved sources with 
i'-z' > 2.2 in the SDSS system are either quasars at z > 5.7 or 
brown dwarfs. The filter transmission profiles and CCD sen- 
sitivity for MegaCam are not identical to those of the SDSS. 
Therefore to determine the optimum colour selection crite- 
ria for high-redshift quasars with CFHT, we have performed 
simulations of the expected colours of quasars in the CFHT 
filters. 

Previous methods to simulate the colours of high-redshift 
quasars have used composite quasar spectra or power-law + 
emission line spectra as the input quasar template. To obtain 
a realistic spread in quasar properties, we use the observed 
spectra of a sample of z ~ 3 quasars to simulate the colours of 
quasars at higher redshift. A similar method has recently been 
presented by Chiu et al. (2005). However, their analysis was 
for the SDSS filters and differed from ours in several aspects, 
so we outline our method in detail here. 

Since our high redshift quasar selection is based on the i',z' 
and J filters, we require optical quasar spectra which sample 
the rest-frame wavelengths of 5.5 < z < 6.7 quasars in these 
filters. This wavelength range can be achieved with the SDSS 
DR3 spectroscopy of quasars at 3.1 < z < 3.2. At these red- 
shifts, there is very little colour incompleteness in the SDSS, 



due to the fact that the quasar colours lie far from the stellar 
locus (Richards et al. 2002). This is important since we are 
using these quasars to assess the efficiency of our own colour 
selection criteria. There are 180 quasars in the DR3 catalogue 
satisfying these criteria. The lack of evolution in quasar spec- 
tra from z = 6 to lower redshift (Fan et al. 2004) justifies our 
use of lower redshift templates. We assume the range of emis- 
sion line strengths, broad absorption line occurrence and spec- 
tral shapes of the z = 3 quasars are representative of the z = 6 
quasar population. 

The only significant difference between the rest-frame UV 
spectra of z = 3 and z = 6 quasars is the much stronger HI ab- 
sorption at higher redshift. Therefore the spectra of the z = 3 
quasars need to have their Lyman forest absorption properties 
changed to be consistent with being located at a higher red- 
shift. The Lyman forest absorption over the redshift range 
2 < z < 6 has recently been studied using high-resolution 
spectroscopy by Songaila (2004). The mean transmission for 
Ly a as a function of redshift is taken from their eqn. 3. This 
relation only applies at z > 4, so the lower redshift transmis- 
sion was determined by fitting to the data tabulated in that pa- 
per. The Ly (3 optical depth was determined using Tp = 0.40r Q 
as found by Songaila. First of all it is necessary to statistically 
correct for the Ly a and Ly (3 absorption already present in the 
spectra of the z = 3 quasars. Ly 7 and higher order absorption 
are not considered since they have a very small effect on the 
flux integrated over any of our filters at any redshift. We also 
set all flux below the quasar Lyman limit to zero. After cor- 
recting for the lower redshift absorption, absorption at higher 
redshifts assuming the quasar is located at 5.5 < z < 6.7 is 
applied. For each z = 3 quasar, we introduce scatter in the 
Lyman forest absorption as a function of redshift (assuming 
the quasar lies at 5.5 < z < 6.7). The amount of scatter in 
Sz = 0.1 bins is estimated as a function of redshift from the 
observations of Songaila (2004). 

To determine the colours of quasars in our filters we con- 
volve the 180 simulated quasar spectra with the i',z' (Mega- 
Cam) and J (CFHT-IR) filters. The filter profiles include the 
effect of CCD quantum efficiency and atmospheric transmis- 
sion. All magnitudes are determined on the AB system. Note 
this is different from the z > 5.7 quasar survey in the SDSS 
where i' and z' are on the AB system, but J is on the Vega sys- 
tem 5 . The colours are determined for a range of redshifts from 
z = 5.5 to z = 6.7. We also ran the simulated quasar spectra 
through the SDSS filter system to determine the SDSS colours 
of our quasar templates. We found good agreement between 
the simulated i' - z' colours as a function of redshift and the 
observed colours of the 12 z > 5.7 SDSS quasars. This pro- 
vides a check on our Lyman forest absorption technique since 
the i'-z' colours are particularly sensitive to this absorption. 

The i'-z' vs z'-J diagram is shown in Fig.^ The colours of 
all 180 simulated quasars are plotted as a function of redshift 
from z = 5.5 at lower-left to z = 6.7 at upper-right. Note the 
tracks of the quasars in colour-colour space follow a similar 
trend to that shown for a simulated quasar in the SDSS system 
in Fan et al. (2003). However, there is a significant difference. 
The SDSS i' and z' filters are not overlapping, whereas there 
is a small, but significant, overlap for the MegaCam filters. 
The effect of this overlap is to decrease the i' - z' colours of 
quasars, particularly in the range 5.7 < z < 6 where the Lya 
line passes through the overlap region. Therefore we need to 
adopt a lower i'-z' selection criteria than the i'-z' > 2.2 used 

5 Jab = /Vega + 0.90 for the CFHT-IR filter. 
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FIG. 1. — i' —z'vsz'-J in the CFHT filters for simulated quasars and 
literature low mass stars. The coloured curves show the tracks in colour- 
colour space of the 180 redshifted z ~ 3 quas ars a s they are moved from 
z = 5.5 at lower-left to z = 6.7 at upper-right (Sec . 12. 11 . The median colours at 
each redshift are calculated to produce the median quasar track (thick black 
line). Intervals of 8z = 0. 1 along this track are marked with circles. The i' — z' 
vsz' —J colours of low mass stars from the literature are s hown as plus signs 
(M dwarfs), triangles (L dwarfs), asterisks (T dwarfs) (Sec . 12. 2^ . The median 
star colours as a function of type are shown with a thick dashed line. The thin 
solid line denotes the high-redshift quasar selection region given in eqns. 1 
and 2. Dotted lines indicate if— z' = 1.5 and f — z! = 1.7. 



by the SDSS if we are to select quasars at similar redshifts. 
The box bounded by the regions 

i'-z' > 1-5 (1) 

i'-z'-l.5(z'-J)>0.6 (2) 

defines the default colour-colour high-redshift quasar selec- 
tion criteria. It can be seen from inspection of the figure that 
these criteria select almost all quasars in the range 5.8 < z < 
6.4. 

To determine the quasar selection efficiency as a function of 
redshift more quantitatively, we have determined the fraction 
of quasars as a function of redshift which fall within this box. 
To account for the fact that the colours of objects detected in 
CCD imaging have a photometric uncertainty we model the 
location of each quasar at each redshift in colour-colour space 
as a 2D gaussian probability density distribution with a = 0.2. 
This error i n col ours is typical for the faintest sources in our 
survey (Sec. l3.2> . We then proceed to calculate the fraction 
of the probability that falls within the selection region. This 
process is repeated for three different i'—zf cuts of 1 .5, 1 .6 and 
1 .7 to determine the effect of imposing different colour cuts 
on the selection efficiency as a function of redshift. 

The results of this process are shown in Fig.|2] Quasars 
can be selected over the redshift range 5.7 < z < 6.5 at 
a completeness fraction of > 0.5. The completeness is 
close to unity at 5.8 < z < 6.4. The effect of imposing a 
stricter cut of i' — z' > 1.6 or 1.7 is fairly small and only has 
an impact at the lower end of the redshift range (up to z = 5.9). 




5.6 5.8 6.0 6.2 6.4 6.6 
Redshift z 



FIG. 2. — The fraction of quasars selected by the colour selection criteria 
of eqns. 1 and 2 as a function of redshift (soli d line). This takes into account 
photometric scatter as described in Sec. 12.11 The dashed and dotted lines 
show the fractions if replacing eqn. 1 withi'— z' > 1.6 and 1.7, respectively. 



2.2. Ultracool stars 

To determine the colours of ultracool stars (M, L and T 
dwarfs) in our filters we convolve observed spectra obtained 
from Sandy Leggett's online database 6 with the i', z' (Mega 
Cam) and J (CFHT-IR) filters. Only stars without significant 
gaps in spectral coverage at the wavelengths of these filters 
were used. The spectra are from the following sources: Tsve- 
tanov et al. (2000), Burgasser et al. (2000), Kirkpatrick et al. 
(2000), Reid et al. (2001), Leggett et al. (2002) and Knapp et 
al. (2004). Fig.[0shows the dwarf colours as plus signs (M), 
triangles (L) and asterisks (T). A sequence is also plotted by 
binning the dwarfs by type and evaluating the median colour 
in each bin. This figure demonstrates that the i' —z' vs z 1 -J 
diagram provides a very clean separation between quasars and 
brown dwarfs for i' — z! >1.3. The i' — z' > 1.5 selection cri- 
teria which we adopt nominally selects dwarfs cooler than 
L2/L3. With significant noise in our colours and the much 
larger population of late-M dwarfs in a flux-limited sample, a 
significant number of slightly warmer stars may be observed 
to have i'—zf > 1.5. 

3. OBSERVATIONS 

3.1. CFHTLS Imaging 

In this paper we use optical imaging from the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) 7 to 
search for high-redshift quasars. The data used come from 
the first official release: T0001. This release includes data in 
the four Deep Survey fields obtained during the period June 
2003 - July 2004. The Deep Survey has a 5-year lifetime, so 
the final dataset will go ~ 1 mag deeper than the data used in 
this paper. 

Full details of this dataset are available elsewhere 8 , but here 
we summarise the details relevant to our quasar search. Each 
of the four Deep fields (named Dl-4) is one dithered point- 
ing with the MegaCam camera. Therefore the total area of 
each field is w 1 sq. deg. The total integration time per field in 

6 ftp://ftp.jach.hawaii.edu/pub/ukirt/skl/ 

7 http://www.cfht.hawaii.edu/Science/CFHTLS 
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the i' and zf filters respectively are Dl: 52.0ks, 12.2ks; D2: 
18.5 ks, 10.1 ks; D3: 59.6 ks, 15. Iks; D4: 58.8 ks, 26.6 ks. 
The seeing FWHM is rj 0.9" in all the images. The astrome- 
try of the images in the two bands are aligned to better than a 
pixel (0.186"). 

To find very red objects in this dataset we generated cata- 
logues of objects detected in the z'-band images using the Sex- 
tractor software (Bertin & Arnouts 1996). We ran Sextractor 
in "double-image" mode to determine /'-band measurements 
for all the objects detected in zf. Magnitudes in both images 
were measured in 2" circular apertures. Aperture corrections 
to total magnitudes were applied for unresolved sources. The 
magnitudes have been corrected for galactic extinction using 
the maps of Schlegel, Finkbeiner & Davis (1998). Note that 
stars may not lie behind the total galactic obscuring column, 
but since the main goal of this paper is quasar selection and 
the applied galactic extinction corrections are smaller than our 
photometric errors, we do not attempt any further corrections 
for stars. 

From these catalogues we drew a sample of objects satisfy- 
ing the colour criterion of eqn. ^ We restricted the selection 
to sources which have z'-band photometric errors er(z') < 0. 10 
to ensure we do not spend time following up spurious detec- 
tions. Although images of these fields in the u*,g' and r 1 fil- 
ters are also part of the T0001 release, we did not impose any 
other colour criteria since the i' images have the longest inte- 
gration time so most selected objects are undetected at shorter 
wavelengths. Every candidate was inspected by eye to ensure 
that it is real. There were many spurious sources, mostly lo- 
cated near to very bright stars, where the diffraction patterns 
in the i' and z' images were different. After excising these 
from our target list we are left with 7, 12, 11 and 5 objects 
satisfying eqn. [Qand ct(z') < 0.10 in fields Dl, D2, D3 and 
D4, respectively. 

3.2. Near-infrared imaging 

As shown for the SDSS by Fan et al. (2001) and dis- 
cussed in Sec.|2] the most efficient way of discriminating be- 
tween high-redshift quasars and brown dwarfs is by using 
near-infrared /-band photometry. Brown dwarfs and very 
low mass stars have red z'— J colours and quasars have blue 
Z 1 —J colours. The exact colour selection criteria used in th e 
CFHQS to separate stars and quasars was given in Sec . 12. II 

Observations at /-band were performed using the CFHT-IR 
camera at the CFHT. CFHT-IR is a HgCdTe 1024 2 array with 
pixel size 0.21 1 ". Observations were carried out in photomet- 
ric conditions on UT 2003 November 9, 10, 2004 May 4 and 
2004 November 23, 24. Our first observations with CFHT-IR 
came a year before the T0001 release of the CFHTLS Deep 
data. For the first two observing runs we used preliminary 
reductions of the CFHTLS dataset generated by one of us 
(SDJG) and also by David Balam and Chris Pritchet. This led 
to the observation of some objects which now have i'-z 1 < 1 .5 
in the most recent (and deepest) reduction. Due to constraints 
imposed by scheduling and poor weather it has not been pos- 
sible to observe a complete sample of i'-z 1 > 1.5 sources to 
a consistent magnitude limit in all 4 fields. We prioritised the 
reddest and brightest sources so that we can quantify our in- 
completeness. The sources not observed all lie close to the 
i'-z' =1.5 boundary and therefore we can account for the 
possibility th ey m ay be quasars at z ~ 5.8 in our selection 
function (Sec.BlV 

8 http://terapix.iap.fr 
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FIG. 3. — Very red objects selected from the CFHTLS Deep i' and z' imag- 
ing which were observed with CFHT-IR at /-band are shown by squares with 
lcr error bars. Filled squares show those objects belonging to the complete 
sample of candidates defined by — z 1 > 1.5, cr(z') < 0.10 and z 1 < 23.35. 
Open squares show sources observed which fail one or more of these criteria. 
Also shown are the median colours of simulated quasars (thick solid line) and 
low mass stars (thick dashed line) from Fig.0 The thin solid line denotes the 
high-redshift quasar selection region given in eqns. 1 and 2. Dotted lines 
indicate i' — z' = 1.5 and ;' — z' = 1.7. 



A total of 24 sources were observed at /-band of which 20 
have i'-z' > 1.5 (the rest have 1.45 < i'-z' < 1.5). 15 out of 
the 35 sources identified in Sec. 13.11 were not observed (five 
in D2, nine in D3 and one in D4). The poor completeness 
in D3 is a consequence of bad weather during the CFHT-IR 
run in May 2004. /-band magnitudes on the AB scale were 
measured using 2" circular apertures and aperture corrections 
applied to give total magnitudes. 

The resulting photometry is given in Table 1 . The i' - z' 
vs zf —J colours of these objects are plotted in Fig.|5] As is 
immediately clear in this diagram, the colours of all 24 tar- 
geted objects are consistent with stars and inconsistent with 
quasars. The three closest to the quasar selection region with 
z' —J < 1 are not part of the complete sample. Two of these 
have i'-z! < 1.5. The other has ct(z') > 0.10 and a large un- 
certainty on its /-band magnitude and hence zf —J colour. It 
is possible that this object is a quasar, but given the poor con- 
straint on its colour and the fact that stars are many times more 
common than quasars at this i'-z' colour, the prior on the like- 
lihood is that it is a star. 

The two reddest objects (at zf - J ~ 3) do lie close to the 
expected colours of z > 6.8 quasars. However such a high 
redshift would require very high-luminosity quasars since at 
this redshift the Lyman-a line is moving out of the z'-band. 
The luminosity of such quasars would be comparable to those 
in the SDSS at z ~ 6 which have a surface density of « 1 
per 500 square degrees (Fan et al. 2004). The probability 
of finding such a quasar in 3.8 square degrees is very low. 
There is a m uch greater probability that they are T dwarfs 
(see Sec. l3.3t . Even if these objects were quasars it would not 
affect the results of the rest of this paper because we only use 
the survey volume for quasars out to z = 6.4. 
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TABLE 1 
Table of observed objects. 



Field 


Name 


i' (AB) 


z' (AB) 


J (AB) 


Complete? 


Dl 


PFT4TT S 1077506 79-047958 4 


25 41-1-0 18 


3n -1- 09 

Lj.JU ZIZ U.Lf^ 


71 64 4- 16 

Z. 1 . UH- ZTZ U. 1U 




Dl 


CFHTT <s 1077456 77-045Q05 Q 


74 30 4-0 10 


97 on 4. n 07 

Zl.OU ZIZ O.v / 


91 514-nin 

Zl.Jl ZIZ V7. i VJ 




Dl 


TFHTT ST077714 12-045723 1 


24 35 + 10 

£rr . J J l \J . 1 \J 


22 81 -t- 06 

ilL.Ol 1 VJ . VJU 


91 40 4- 08 

1 .^7 1 VJ . VJ O 


/ 


Dl 


TFHTT S T0774SS 66-040607 7 


91 70-Lnm 
z. i . / o r_r_ vj . vj j 


70 17 _|_ n 

Z.U. 1 / ZLL VJ . VJ J 


IS 1 ft 4- 04 




Dl 


TFHTT ^ T07741 S SQ 0404T1 1 


7a SO 4- 07 


77 71 _|_ 04 
Z.Z..Z.1 in u.u i + 


70 KS 4- 07 


V 


Dl 


TFHTT 's T077SSQ 1Q-040400 


70 47 4_ n nft 
_r_ vj.l/u 


71 04. _i_ 01 


70 11 4- 04 

Z.VJ. 1 1 IC VJ . U't 




Dl 


TFHTT S T07777S OQ-04011 S 7 


7 a 40 4- 06 


71 Q£ _L 04 
Z1.7U IC VJ . yji 


70 44 4- 04 




D2 


CFHTT S 11001 13 06+022622 4 

v I 111 l_j*J J 1 \JV7 1 1 _J .\JVJ~VJZ*Z.\JZ.Z*.*"T 


> 26.35 


22 70 4. n 07 


1 q 71 4- 04 
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D2 


CFHTLS JT 00132.04+020022.4 


22.91 ±0.05 


21.39 ±0.03 


19.73 ±0.05 


/ 


D3 


CFHTLS J 1 4 1 659.9 1 +52 152 1 .9 


23.89 ±0.04 


22.32 ±0.03 


20.07 ±0.08 


■/ 


D3 


CFHTLS J 1 4 173 1 .77+523920.4 


23.73 ±0.04 


22.17 ±0.03 


20.13 ±0.07 




D4 


CFHTLS J221705.05-172932. 5 


24.47 ±0.10 


22.67 ±0.05 


20.77 ±0.11 




Dl 


CFHTLS J022757.64-043019. 5 


25.03 ±0.16 


23.35 ±0.12 


22.53 ±0.30 


X 


Dl 


CFHTLS J0225 16.57-041402.9 


24.87 ±0.14 


23.38±0.11 


22.19±0.35 


X 


Dl 


CFHTLS J0226 14.42-044600.4 


24.26 ±0.10 


22.81 ±0.08 


21.91 ±0.32 


X 


Dl 


CFHTLS J022734.82-042438.7 


23.52 ±0.06 


22.05 ± 0.04 


21.10±0.16 


X 


D2 


CFHTLS J0959 1 4.80+023655 .2 


> 26.35 


23.41 ±0.10 


20.45 ± 0.07 


X 


D4 


CFHTLS J221430.53- 1 80230.5 


25.11 ±0.16 


23.41 ±0.09 


> 21.98 


X 


D4 


CFHTLS J22 1 457.74- 1 72 1 04.2 


25.32 ±0.17 


23.66 ±0.10 


> 21.98 


X 


D4 


CFHTLS J221437.22-1 80417.4 


24.08 ±0.08 


22.59 ±0.05 


21.03 ±0.09 


X 



NOTE. — Objects selected by their red i — z colours from the CFHTLS Deep fields which were observed at /-band with CFHT-IR. The magnitudes listed are total magnitudes 
derived from aperture magnitudes. The final column indicates whether the object satisfies the i and z selection criteria which form the basis for our complete sample of quasar 
candidates for /-band imaging. These criteri a are i' -z > 1.5, a(z') < 0. 10 and z < 23.35. The three-colour photometry indicates that all these objects are ultracool dwarfs fFigj3l. 
Their properties are discussed further in Sec J3.3l 



3.3. Ultracool dwarfs detected 

The two reddest objects (CFHTLS J095914.80+023655.2, 
i'-z' > 2.9 and CFHTLS J 1001 13.06+022622.4, i'-z' > 
3.5) both have z'-/~3 and are T dwarfs. The constraint 
from the i'-z' and z' -J colours indicate that the former is 
cooler than Tl and the latter cooler than T4. To charac- 
terise their type more effectively we obtained //-band pho- 
tometry of them with CFHT-IR. Both objects are very blue 
at J-H confirming their T dwarf nature (even bluer than 
a z - 7 quasar). CFHTLS J095914.80+023655.2 has J- 
H(AB) = -0.26 ±0.09 and CFHTLS J1001 13.06+022622.4 
has J-H(AB) = -0.70 ±0.08. Using the transformation 
J - //(Vega) = J - //(AB) + 0.44 for comparison with the 
plot of J-H vs spectral type in Leggett et al. (2002), 
we find that CFHTLS J095914.80+023655.2 is type T3 and 
CFHTLS J1001 13.06+022622.4 is type T6. With J AB « 20, 
they are 3 -4 magnitudes fainter than most SDSS and 2MASS 
T dwarfs. Approximate distances based upon the relationship 
between absolute magnitude and spectral type (Tinney, Bur- 
gasser and Kirkpatrick 2003) are 150 and 70 pc, respectively, 
making them two of the most distant T dwarfs known. 

Only four or five of the analysed objects have colours typi- 
cal of L dwarfs (i'-z' > 1 .5 and z' -J - 2.0). The z'-J < 1 .7 
colours of the rest indicate that they are late-M dwarfs, with 
a true i'-z' below 1.5, scattered into our i'-z' box by 1-3 a 
random noise excursions (this scattering is also apparent in the 
SDSS sample, see e.g., Fig. 3 of Fan et al. 2003). Within the 
limited statistics, the number of detected L dwarfs is consis- 
tent with an extrapolation of the densities seen in the shallow 
whole sky surveys. 

3.4. Completeness 

To turn the non-detection of quasars in this dataset into a 
useful limit on the quasar space-density and luminosity func- 



tion, it is necessary to evaluate the completeness of our survey. 
There are two relevant selection criteria which will affect the 
completeness as a function of redshift and absolute magni- 
tude: 

• The z'-band magnitude limit as a function of the survey 
area 

• The colour selection criteria 

Excluding the field edges which we have masked out of our 
catalogues, the z'-band images show a fairly uniform noise 
and therefore detection limit across the field. Since the major- 
ity of our candidates have z' < 23.35 and point-sources of this 
magnitude have cr(z') < 0.10 across most of the survey, we 
adopt a magnitude limit for the present analysis of z' < 23.35. 
To determine the efficiency with which objects could be found 
as a function of magnitude we ran simulations by inserting 
10000 artificial stars into each z'-band image and using Sex- 
tractor to recover them. Even at bright magnitudes, about 5% 
of sources are not recovered because they overlap with other 
objects or regions affected by bright stars. We note that this 
fraction is much lower than the w 20% of the area in the mask 
region file released by Terapix. We find their masking to be 
much too conservative for our purpose. By z' = 23.5 the re- 
covered fraction has decreased to « 80% in all the images. 
For each field we determine the area in which sources can be 
found as a function of their magnitude. At bright magnitudes 
the total area surveyed is 3.83 sq deg, dropping to 3.32 sq deg. 
at our limit of z' = 23.35. 

The colour selection criteria were discusse d ex tensively 
with simulated quasar and stellar colours in Sec. 12. II Here we 
recap these selection criteria and evaluate the extent to which 
colour selection impacts our completeness. The nominal 
colour selection criteria are i'-z' > 1.5 and i'-z' -1 .5(z'-J) > 
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0.6. No observed objects lie within this box. Since we were 
unable to observe all the candidates with i'—zf > 1.5 down 
to the magnitude limit, our completeness level in terms of the 
i' —z' colour varies from field-to-field. Field D 1 is complete to 
i'-z' > 1.5, D2 and D3 are complete to i'-z' > 1.7 and D4 is 
complete to i'-z' > 1.6. Note that even though only two out 
of eleven candidates in D3 were observed, all the unobserved 
objects have 1.5 < i'-z' < 1.7. 

Fig-El s howed the completeness due to colour cuts as deter- 
mined from the simulated high-redshift quasar colours. The 
effect of implementing redder cuts than the nominal i'-z' > 
1 .5 is to exclude quasars in the range 5.7 < z < 5.9. The com- 
pleteness at z = 5.8 is 89%, 80% and 68% for i'-z' > 1.5, 1.6 
and 1.7, respectively. At z = 5.9, the colour completeness 
for all three cuts are > 98%. Therefore using a redder cut 
in i'-z' makes only a very small difference to our sensitivity 
to quasars across the full range of 5.8 < z < 6.4. This means 
that we can include fields like D2 and D3 which have not been 
so well followed up in the near-infrared whilst losing only a 
small amount in the completeness of high-redshift quasar se- 
lection. 

The magnitude and colour completenesses are combined to 
generate an effective volume of our high-redshift quasar sur- 
vey as a function of redshift (over the interval 5.8 < z < 6.4) 
and absolute magnitude. Absolute magnitudes of the quasar 
continuum at a rest-frame wavelength of 1450A are calcu- 
lated for the possible range of apparent z'-band magnitudes. 
These calculations use ^-corrections generated from the mean 
^-c orrec tions of the 180 redshifted z ~ 3 quasars discussed in 
SecO 



4. CONSTRAINING THE QUASAR LUMINOSITY FUNCTION 

We now use the non-detection of high-redshift quasars in 
the CFHTLS Deep fields to constrain the luminosity function 
(LF). We parametrise the LF using the double power-law form 
that provides a good fit at lower redshift (Croom et al. 2004), 
i.e. 



$(Mi450,z) = 



(3) 



lQQA(a+l)(M 14 50-M^ 50 ) + iQOA(l3+\)(M 145 o-M^ 5B ) ' 

The luminosity distribution of SDSS z ~ 6 quasars from Fan et 
al. (2004) is used to fit the bright-end normalisation $(Mi45o = 
-26) as a function of the bright-end slope j3. Taking the range 
in (3 allowed by the SDSS and the value of $(M I4 5 = -26) at 
each /?, we set up an array of possible double-power-law LFs. 
The faint-end slope a will not be strongly constrained by our 
data and therefore we do not treat this as a free parameter. 
Instead we assume that a takes a value of -1.0,-1.5 or -2.0 
and repeat our analysis for each of these values. The best-fit 
values at the highest redshifts at which a has been measured 
(z ~ 2-3) lie in the range -1.1 to -1.8 (Boyle et al. 2000; 
Croom et al. 2004; Hunt et al. 2004; Richards et al. 2005). 
We will consider the effect of a later on. Sec. 13.41 discussed 
our survey volume as a function of redshift and absolute mag- 
nitude. We use the array of possible LFs combined with 
the survey volume to calculate the number of 5.8 < z < 6.4 
quasars /i that we would expect to have discovered in our sur- 
vey as a function of (3 and break magnitude M* 450 . The prob- 
ability of not finding any quasars given this LF is then simply 
given by the Poisson probability P = e~^. 

The result of this analysis is shown in Fig.|4] Curves are 
plotted at confidence levels of 0.32,0.05 and 0.01. Sep- 
arate panels show the different results obtained by assum- 
ing faint-end slopes of a = -1,-1.5,-2. Since our survey is 
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FIG. 4. — Constraints on the z = 6 quasar luminosity function from a va- 
riety of optical surveys. Each panel shows the constraints on the bright-end 
slope /3 and break magnitude M* 450 . The three panels assume values for the 
faint-end slope of a = -1 (upper), —1.5 (middle) and -2 (lower). Confidence 
intervals are drawn for the CFHQS, SDSS and the CDF-N (Sec.|4). The best 
fit parameter values for the 2dF at z = 2 (Croom et al. 2004) are shown for ref- 
erence. The grey shaded region shows the values of parameters necessary for 
the z = 6 quasar population to produce enough ionizing photons to reionize 
the universe as described in SeclSI 



only sensitive to quasars brighter than M1450 = -23, our con- 
straints depend weakly on a unless the break magnitude is 
very bright (M'^ 45[) < -25). The evolution in the 2dF quasar 
LF is fit by pure-luminosity evolution with the break evolving 
from Mf 4sn = -21 at z = to Mf 4V , = -25 at z = 2 (Croom et 
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FIG. 5. — Plausible double power-law parametrisations of the z = 6 quasar 
optical luminosity function are shown as colored lines. Only those models 
with confidence levels > 0.05 based on the CFHQS and SDSS constraints are 
shown. These models are normalized to the SDSS results of Fan et al. (2004) 
which is why they have approximately the same density at M1450 = -27. The 
luminosity function at z = 2 from the 2SLAQ + 2dF surveys (Richards et al. 
2005) is shown as a black line. Our work finds that the large decrease in the 
space density of luminous quasars from z = 2 to z = 6 found by the SDSS also 
extends to lower luminosities (M^o ~ -24). 



al. 2004). Since the number density of luminous optically- 
selected quasars declines at redshifts beyond 2-3 (Fan et 
al. 2004), continued pure-luminosity evolution would lead 
to a fainter M* 450 at z = 6 than at z = 2. Meiksin (2005) de- 
termined that a model with continued pure-luminosity evo- 
lution that fits the luminous SDSS z = 6 sample would have 
M* 450 = -22.9 at z = 6.0. Although we do not necessarily ex- 
pect the high-redshift evolution to be characterised by lumi- 
nosity evolution, it is unlikely that the break magnitude would 
be brighter than it is at z = 2, so we consider M* 450 = -25 as a 
likely lower limit. 

Barger et al. (2003) reported a negative search for z > 5.2 
quasars in the Chandra Deep Field-North (CDF-N). Their 
work goes deeper than our study (to z' = 25.2) but covers 
a much smaller area (0.031 sq. deg). We have therefore 
repeated our analysis of LF constraints using the lack of 
5.8 < z < 6.4 quasars in the CDF-N. The 0.05 confidence con- 
tour on (3 and M\ A5Q from the CDF-N is plotted on Fig.|4] It 
can be seen that our constraint from CFHQS is considerably 
stronger than that from the CDF-N due to our survey volume 
being w 100 times greater. Also plotted on Fig.[4]are the con- 
straints on (3 from the luminosity distribution of the SDSS 
(Fan et al. 2004). The CFHQS provides much better lower 
limits on (3 than the SDSS. For all the values of a consid- 
ered, the 0.05 confidence limit is > -3.2 if M* 450 = -24 and 
13 > -3.1 ifM* 450 = -23. 

FigHJshows the LF constraints from the CFHQS more di- 
rectly. The colored curves correspond to LF models which 
are allowed at > 0.05 confidence. We can use the constrained 
LF parameters to place an upper limit on the space-density 
of quasars at z = 6 brighter than a certain absolute magni- 
tude limit. We adopt a limit of Mi 4 5o < -23.5 since at this 




-1.5 -2.0 -2.5 -3.0 -3.5 -4.0 -4.5 



FIG. 6. — The ionizing photon production rate density for the z = 6 quasar 
population as a function of /3 and M* 4SQ . Contours are drawn for three values 
of the luminosity function faint-end slope: a = -1 (solid), —1.5 (dotted) and 
-2 (dashed). Contours are plotted at 10 49 , 10 50 , 10 51 , 10 52 photons s" 1 Mpc~ 3 . 



magnitude cut there is almost no dependence of the den- 
sity upper limit on a,/3, or M* 450 . We find that p(Mj 4 go < 
-23.5) < 1.7 x 10 _7 Mpc" 3 at 95% confidence. This limit is 
a factor of 25 lower than the density at z = 2 measured in the 
2SLAQ + 2dF surveys (Richards et al. 2005). This confirms 
that the decline in the space-density of luminous quasars at 
z ~ 6 (Fan et al. 2004) is mirrored by the lower luminosity 
population. A decline in the space-density of low luminos- 
ity quasars has been seen in the range z ~ 3-5 (Barger et al. 
2003; Wolf et al. 2003; Cristiani et al. 2004; Sharp et al. 
2004), but the large volume and faint magnitude limit of our 
survey extends this result to z ~ 6. 

5. A LIMIT ON THE QUASAR CONTRIBUTION TO THE IONIZING 
BACKGROUND AT Z ~ 6 

If the high-redshift quasar LF can be fit by a double power- 
law with a break at M* 450 ^ - 23 and a > -2, then a large 
fraction of the luminosity emitted by the quasar population 
will be emitted by quasars that we could have found in our 
survey. This means that our survey is useful for providing an 
upper limit to the number density of ionizing photons emitted 
by quasars at z = 6. 

To estimate the number of hydrogen-ionizing photons emit- 
ted by a quasar with a given Mi 4 5o, we assume the typical 
spectral shape determined for a large sample of quasars by 
Telfer et al. (2002). We assume that 100% of these photons 
escape from the host galaxy. The ionizing photon emission 
density is calculated by multiplying the photon emission rate 
per quasar by the number density of quasars and integrating 
over the range of quasar luminosities. We perform the inte- 
gration down to quasars as faint as Mj 4 5o = -18. 

Fig.|5]shows the ionizing photon production rate density for 
the z = 6 quasar population as a function of (3 and M\ A5Q . Con- 
tours are plotted for the three different values of a consid- 
ered in the previous section. Note that the shape of the con- 
tours is quite similar to the shape of the confidence limits on 
these parameters determined from our survey for high-redshift 
quasars. This is due to the reason outlined previously that the 
main constraint from our survey is on the number density of 
quasars brighter than Mj 4 5o = -23 and such quasars dominate 
the ionizing photon budget unless a is very steep. We note 
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that the ionizing photon production rate density in Fig. [6] is 
consistent with the calculations of Fan et al. (2001), but not 
with those of Yan & Windhorst (2004) who found a factor of 
10 lower rate for the same quasar LF parameters. The reason 
for this discrepancy is unclear. 

There is considerable debate about the photon density re- 
quired to reionize the diffuse neutral hydrogen in the universe 
at z = 6. The main uncertainty is due to recombination of ion- 
ized hydrogen atoms, which depends upon the dumpiness of 
the IGM. A detailed discussion of the dumpiness issue is be- 
yond the scope of this paper. We follow Meiksin (2005) and 
assume that the dumpiness leads to 5 ionizing photons per 
hydrogen atom necessary for reionization at z = 6. Note that 
this is substantially fewer than the 30 ionizing photons per 
hydrogen atom used by Madau et al. (1999) which has been 
widely used by other authors. Under our assumptions the ion- 
izing photon production rate density necessary for reionizing 
the universe at z = 6 is 3 x 10 5() photons s" 1 Mpc" 3 . 

The region of /3-Mj* 450 space where the quasar population 
alone can provide the necessary ionizing flux is shown by the 
grey shaded regions in Fig.|4] For any value of a the shaded 
region does not overlap with the CFHQS 0.05 confidence con- 
tours and quasars do not provide enough ionizing photons to 
reionize the universe. For a < -2, the quasar emission is 
dominated by low luminosity quasars and a population with 
a break at M* 450 > -23 could still reionize the universe, al- 
though the details then depend critically upon how far down 
the LF one integrates. All the evidence at intermediate red- 
shifts points towards a > —1.8 (Boyle et al. 2000; Croom et 
al. 2004; Hunt et al. 2004; Barger et al. 2005; Richards et 
al. 2005; but see Hao et al. 2005 for evidence that a ks -2 in 
the local universe) making it unlikely that a could be steeper 
than -2 at z = 6. We conclude that a 95% confidence upper 
limit on the ionizing photon rate density provided by quasars 
at z = 6 is 3 x 10 50 photons s" 1 Mpc" 3 . 

6. CONCLUSIONS 

We have used 3.83 sq. deg. of optical imaging with a mag- 
nitude limit of z' = 23.35 to search for high-redshift quasars. 
Near-infrared follow-up of candidates selected by their red 
i' - z' colours all have z' —J colours consistent with being low 
mass stars: two are T dwarfs, four or five are L dwarfs, and 
the rest are late-M dwarfs scattered into our selection box by 
photometric errors. 

The lack of quasars at z > 5.8 in these data has been used to 
provide new constraints on the luminosity function. Several 
surveys underway at CFHT will image between 100- 1000 sq. 



deg. in multiple filters to somewhat shallower depths than the 
data in this paper. Assuming the bright-end slope is not much 
flatter than the value at z = 4 of (3 = -2.6 (Fan et al. 2001), 
these surveys will allow the identification of many z > 5.8 
quasars giving a much more accurate determination of the lu- 
minosity function and provide targets for studying the ioniza- 
tion state of the IGM. 

We have considered the ionizing photon output at z = 6 and 
concluded that quasars provide < 3 x 10 50 photons s" 1 Mpc" 3 
at 95% confidence and assuming that the faint end slope 
a > -2. The ionizing photon output from the z ~ 6 star- 
forming galaxy population is subject to significant uncer- 
tainties as outlined earlier, but the best estimates give a rate 
> 10 51 photons s" 1 Mpc" 3 (Yan & Windhorst 2004; Bunker et 
al 2004; Stiavelli et al. 2004). Therefore quasars provide 
less than 30% as many ionizing photons as galaxies at z ~ 6. 
The evolution of the neutral hydrogen fraction inferred from 
quasar absorption spectra shows that the universe had been al- 
most completely reionized by z = 6 (Fan et al. 2004; Songaila 
2004). Since something has to keep the universe reionized at 
this redshift, the most straightforward explanation is that the 
galaxy population is responsible. 
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